Vortex matter dynamics in YBa 2 Cu 3 07_ a with BaZrC>3 nanocolumns as revealed by dc 
and high-frequency experiments: role of anisotropy and directional pinning. 
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Measurements of anisotropic transport properties (dc and high-frequency regime) of driven vortex 
matter in YBa2Cu307_a; with elongated strong-pinning sites (c-axis aligned, self-assembled BaZrC>3 
nanorods) are used to demonstrate that the effective-mass angular scaling takes place only in intrinsic 
physical quantities (flux-flow resistivity) , and not in pinning-related Labusch parameter and critical 
currents. Comparison of the dynamics at different time scales shows evidence for a transition of the 
vortex matter toward a Mott phase, driven by the presence of nanorods. 

PACS numbers: 74.25.N-, 74.25.Sv, 74.25.nn, 74.25.Wx 



The role of nanoinclusions in the vortex pinning prop- 
erties of YBa2Cu307-a; (YBCO) films is a matter of 
strong interest 1-4 in view of the great potential of YBCO 
films and tapes with strong pinning (e.g., in coated con- 
ductors). In this field, BaZrC>3 (BZO) as a component 
of nanoinclusions received a particular attention. 1-8 In- 
terestingly, it has been observed that in films grown by 
Pulsed Laser Deposition (PLD), BZO self-assembles in 
the shape of nanorods 1 of typical size ~ 5 nm in diam- 
eter and 30 — 150 nm in length, oriented approximately 
perpendicular to the film plane, that is along the c-axis. 
By contrast, chemical methods tend to yield nanoparti- 
cles instead", of typical size of ~ 15 nm. Such nanos- 
tructures, of typical dimension of the vortex core, are 
thus ideal candidates for strong core pinning. In fact, 
pinning-related properties (critical current density J c , 
pinning force F p = J C B, pinning constant or Labusch 
parameter k p ) are largely enhanced by addition of BZO 
nanoinclusions. l ~' Despite the significant technological 
results, the fundamental vortex physics of anisotropic su- 
perconductors with nanorods is still under debate. 

Aim of this paper is to study the competitive effects 
in the vortex dynamics between the structural, unavoid- 
able anisotropy of the YBCO matrix and the preferred 
orientation introduced by BZO nanorods. By exploit- 
ing the features and the properties of the dynamics at 
different time scales we show that the anisotropic-mass 
angular scaling exists and is limited to intrinsic proper- 
ties, such as the flux-flow resistivity pg. By contrast, 
pinning mechanisms due to the layered structure and to 
the nanorods dominate, so that in pinning-related quanti- 
ties the angular scaling is completely destroyed. Finally, 
we find evidence for a Mott-insulator-like behavior of the 
vortex matter when the magnetic field is aligned with the 
nanorods. This behavior persists with a field misalign- 
ment up to ~ 30° with respect to the nanorods direction. 

A set of YBCO films, 120 nm thick, c-axis oriented, 
were grown on SrTiOs substrates under identical condi- 
tions by PLD from targets with BZO powders at 5% mol. 
as extensively reported elsewhere. 1 ' Transverse TEM im- 
ages showed columnar-like defects, approximately per- 



pendicular to the film plane. 9 Such defects were absent 
in pure YBCO films. The density of columns is consis- 
tent with an equivalent matching field of about 0.85 T. 
T c ~ 90 K (zero dc resistance criterion) was consistently 
found in all samples, with zero-field J c ~ 3.7 MA cm' 2 
at T = 77 K. 

Very different vortex dynamics were studied by the 
dc transport critical current density J C {H,9), and the 
high-frequency (48 GHz) ac transport measurements of 
the field-increase of the complex resistivity Ap(H, 0) = 
p(H, 9) — p(H = 0) = Apt (H, 9) + iAp 2 (H, 8) . Here, 9 is 
the angle between the applied field H and the c-axis. 

Microwave measurements were performed on unpat- 
tcrncd samples I and II by a sapphire cylindrical dielec- 
tric resonator operating at ~48 GHz 1 " and modified to 
work in transmission. Special care was taken 11 in order 
to avoid measurements in correspondence to the well- 
known substrate resonances arising in SrTiOs. Field 
sweeps (p,oH <0.8 T) were performed at different angles, 
as well as angular rotations at fixed fields. Due to the 
small signal when — > 90°, we performed measurements 
at T ~ 80 K, sufficiently below T c to avoid contributions 
to the microwave response from pairbreaking, but still 
with an acceptable signal-to- noise ratio. Fig. la reports a 
typical field sweep Ap(H, 9 = 0°), Fig.s lb and lc report 
angle-rotations at poHi = 0.4 T and poH 2 = 0.6 T. 

Measurements of the transport critical current density 
J C (H,9) were taken on sample II after patterning as a 
strip 30 /im wide and 1 mm long, in the standard four- 
contact configuration, with a 1 pN cm -1 criterion. 

In YBCO with BZO nanorods we expect three sources 
of anisotropy in the vortex dynamics: (a) the mass 
anisotropy leads to larger upper critical field H c2 as H is 
tilted away from the c-axis. Thus, at fixed field we ex- 
pect a reduction of the dissipation solely due to smaller 
reduced field H/H c2 (0). No pinning is conceptually in- 
volved in this effect, (b) the naturally arising intrinsic 
pinning due to the layered structure itself of YBCO as 
a pinning potential, and (c) the strong pinning by BZO 
nanorods, the topic of major interest here. 

The unambiguous identification of the directionality 
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of pinning requires a reliable determination of the elec- 
tronic mass anisotropy involved in (a). This is usually 
done by exploiting the so-called scaling property. 12 ' 13 
Let H C 2(8) = if c2 (O°)/e(0) be the angle- dependent up- 
per critical field, where e 2 (8) 



7 ~ 2 sin 2 I 



cos 



continuous (3D) anisotropic superconductors, and 7 = 
H c2 ($0°)/H c2 {Q o ) is the anisotropy ratio (7 ~ 5 8 in 
YBCO). In the London approximation a thermodynamic 
or intrinsic transport property q depends on the applied 
field H and angle 8 as q(H,6) = s q q[He(0)], where s q 
is a scale factor typically equal to 1, 7 _1 or e ±1 (6') de- 
pending on the observable to be scaled (e.g., s p = 1 for 
in-plane resistivity) . The typical procedure is to look for 
an experimentally determined angular function f(9) such 
that the data for all angles and fields collapse over a sin- 
gle curve when plotted against the reduced field H/ f(8), 
and then compare / to the theoretical expectation. 14 

We stress that this scaling rule is not theoretically 
grounded when sources of anisotropy other than the 
electronic mass affect the observable under investigation 
(e.g.: point defects are likely to leave the scaling prop- 
erties untouched, differently from extended defects). In 
fact, we found that no angular scaling was possible nei- 
ther for Ap(H,8) nor for J c (H,8). ir ' This is not surpris- 
ing: J c is inherently not an intrinsic quantity, and Ap in- 
cludes contributions from free motion of vortices (which 
is an intrinsic phenomenon) and from pinning (extrin- 
sic) . We gain more information by extracting the genuine 
flux flow resistivity pg(H 1 8) and the Labusch parameter 
k p (pinning constant) from the measured Ap, using the 
procedure described with detail previously. 16 This is a 
powerful feature of high-frequency measurements: both 
the flux-flow resistivity and the pinning constant can be 
obtained from the measurements. Since k p = V 2 U eval- 
uated at the pin position (U is the pinning potential), 
it is a measure of the steepness (curvature) of the po- 
tential wells at the bottom. The angular dependence of 
Pff(Hi,8) (i = 1,2) in sample II is reported in Fig. 2a. 
Fig. 2b shows the scaling of the same data over the curve 
taken at 8 = 0°, when plotted as a function of the ap- 
plied field scaled by an experimentally-determined scal- 
ing function f(8) (we discuss f(8) later). For complete- 
ness, Fig. 2d reports the similarly scaled field-sweeps of 
Pff in sample I, performed at several 8. The log scale 
emphasizes the scaling at low fields. It is immediately 
apparent that a nearly perfect scaling is obtained, both 
with field-sweeps and with angle-rotations. 

The experimentally determined angular scaling func- 
tions for field-sweeps, Hi rotations and H 2 rotations are 
reported in Fig. 2c, and clearly describe a unique curve 
that thus does not depend on the sample and on the 
applied field, consistently with an intrinsic phenomenon. 
Thus, as a first result of this paper, we can state that the 
genuine pg obeys the scaling rule, as opposed to J c and 
to the measured, raw Ap. 

In order to compare f(8) with theoretical expressions 
one must take into account that the Lorentz force changes 
with the tilt angle of the magnetic field, due to the cir- 



cular microwave currents in the experimental setup. 1 " 
Working within the 3D anisotropic Ginzburg-Landau 
theory, taking p c ,n/pab,n — 7 2 and averaging over the 
circular current pattern, and taking into account the 
experimental feature of our data, that pff(H) oc H@ 
with /3 ~ 0.8, we find that our measured flux flow 
resistivity is still expected to fulfill the scaling rule: 1!l 
Pff(H,8) — pff(H/f(8)), with the expected scaling func- 
tion f{8) ^ e -1 (#) due to the varying Lorentz force con- 
tribution. We get: 19 
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where the term in square brackets is the correction due to 
the Lorentz force. Eq.(l) has no fit parameters, since the 
experimental value at 8 = 90° fixes 7 and, as shown in 
Fig. 2c, it is in excellent agreement with the data. We get 
7 = 5, within the range of accepted values for YBCO. 2 " 
We stress that the effects of the mass anisotropy and of 
the Lorentz force on the anisotropic response have been 
obtained from the analysis of pg alone, and they are not 
subjected to further adjustments. 

We now turn to the issue of the directionality of pin- 
ning, using the data taken with increased angular accu- 
racy in sample II. To obtain the angular dependence of 
k p (H,8) from the data, one has to remove the Lorentz 
force angular contribution, which enters in the measured 
quantity analogously to what happens for pg . Since from 
Eq.(l) one can explicitly extract such correction, the true 
k p (H, 8) can be obtained from the data 16 ' 19 . 

In Fig. 3b we plot k p (Hi,8) (i = 1,2): a large peak 
can be seen at 8 = 90°, where both the mass anisotropy 
and the a — b plane pinning are expected to decrease the 
dissipation. No directional effect acting around 8 ~ 0° 
seems evident. These observations do not imply that the 
anisotropy of k p (8) is dictated by the mass anisotropy: 
for this to be true, the so-called BGL scaling 12 should 
apply, and s^(9)k p (H,8) vs. He(8) would describe a 

single curve. Since s^ 1 (8) = e(8), 12,19 we plot in Fig. 3a 
the data of Fig. 3b as k p (H \8)e(8) vs. He(8), together 
with k p (H, 0°) as derived from the field sweep with H || c: 
the scaling clearly fails, the main origin of the anisotropy 
of kp is not the mass anisotropy. The angular depen- 
dence in Fig. 3a necessarily originates from other effects, 
demonstrating the existence of directional pinning unre- 
lated to the anisotropic mass, stronger when the field is 
aligned with the nanocolumns (large He(8) in Fig. 3a) or 
with the a-b planes (8 — > 90° and small e(8) in Fig. 3a). 

The comparison of the pinning constant, as obtained 
from microwave measurements, and the dc critical cur- 
rent density J c , elucidates further the physics of vortex 
matter in YBCO with nanocolumns. The maximum pin- 
ning force per unit length is k p r p = J c $o where r p is the 
pinning range of a defect, and by very general arguments 
(e.g., Ref. 21), r v ~ £ for core pinning. In inclined fields 
one has 22 ' 23 r p (8) ~ £ab e (8), where £ a b is the a-b plane 
coherence length. Thus, from microwave data we define 
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an equivalent critical current density as: 



c. short 



(if, i 



k p {H,e)u<e) 



(2) 



where c ~ 1, and the subscript "short" indicates the very 
important point that the dc and microwave vortex dy- 
namics differ substantially: J c is measured in the regime 
where ffux lines are depinned by a sufficiently strong di- 
rect current, and steady motion arises. J c ,short is defined 
in the subcritical microwave currents regime, where the 
microwave field induces only very short (~ 1 A) oscilla- 
tions of the vortices around their equilibrium position. 

Comparison of J c , s hort to tric actually measured dc 
J c is shown in Fig. 3c. Common features are a maxi- 
mum centered at 0°, a local minimum around 85° and 
a steep peak at 90°. The latter clearly originates from 
a-b planes pinning whereas the former is ascribed to the 
correlated pinning due to BZO inclusions. 24 It should be 
noted that the maximum at 0° is extremely wide, thus 
suggesting that BZO nanorods are strong pinning cen- 
ters even in tilted magnetic field. Moreover, with reason- 
able U{T) = Wv 71 ~ (T/T c ) 2 , where £ = 12 A, and 
c ~ 0.16, the short-range-dynamic (J c . short) and long- 
range dynamic ( J c ) curves exactly coincide in the inter- 
mediate angular range 35° < 9 < 80°. Thus, in the 
angular range where flux lines are most likely segmented 
between nanorods, the angular dependences of J c and 
Jc. short arc identical, despite the very different dynamics. 
However, J c ,short{9) and J c (9) strongly depart from each 
other as the field aligns with the nanorods: for 9 < 35° 
Jc. short shows only a broad hump, over which a second, 
steeper hump appears in J c . The further enhancement 
of J c (9) is then a dynamic effect: the largest peak at 
= 0° appears only when large vortex displacements are 
involved (dc). Thus, we can state the important result 
that nanocolumns enhance J c by means of two distinct 
mechanisms: strong core pinning, which gives rise to the 
broad hump in J c (9), and the reduced vortex mobility ef- 
fect, which gives rise to the additional peak in J c (9) for 
< 35°. The latter effect is a manifestation of the pre- 
dicted Mott-insulator phase for fluxons, 25 as we argue in 
the following. Essentially, each vortex is pinned by an 
extended defect and when it is forced to move away from 
it by the external current, its mobility is drastically re- 
duced by the lack of free sites where to move. When the 
field is tilted away from the nanocolumns, fiuxons are 
no longer strongly pinned by the extended linear defects 
for the whole length, so that the Mott-insulator effect 
disappears at a certain angle (~ 35° in our case). De- 
formation of the flux lines (steplike) may yield stronger 
pinning than with point pins only, hence the same wide 
hump in J c .short{9) and J c (9). 

This scenario is fully consistent with the results of 



DC characterizations:" there, we observed a change of 
regime in the pinning force F p approaching the equiva- 
lent matching field (the field where each linear defect ac- 
commodates one flux line). Above that field, the Mott- 
insulator effect vanishes and the pinned vortex matter 
behaves similarly to a more or less disordered lattice in 
presence of strong point pins. Strong extended pinning 
takes place below the pseudo-matching field, as revealed 
by careful analysis of the magnetization relaxation. 26 

We now estimate quantitatively the dynamic effect on 
pinning. A natural way to define the angular dependence 
of the dynamic effect is to introduce the parameter: 



Jc(9) - Jc, short (0) 



Jc 
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When no dynamic effect is present, that is the critical 
current density does not depend on the vortex displace- 
ment, a /j, = 0. From the argument that led to Eq.(2) 
one readily see that can be interpreted as the per- 
centage increase of r p due to the dynamic effect. In Fig. 4 
we plot a ^{9). As it can be seen, a M ^ when the field 
aligns with the directional pinning sources. In particular, 
drops with increasing 9, and vanishes at intermediate 
angles, before a further increase close to the a-b planes. 
The inset of Fig. 4 shows that a M drops exponentially with 
cos 9, the projection of the field along the nanocolumns. 
We found that the following empirical expression 2 ' de- 
scribes very well the angular behavior of 



a^{9) = a^O) 



e Ccos0 _ 2 



(4) 



where £ ~ 14.5 is obtained from the exponential behavior 
for cos 9 > 0.9. The physical significance of the empirical 
fit is that the dynamic Mott-insulator pinning effect is 
a fast vanishing function of the fraction of the flux line 
that can be pinned and "caged" by the extended defect. 

In conclusion, by comparing the angle and field de- 
pendence of vortex dynamics over different time scales 
(dc critical current and microwaves) we have shown that 
angular scaling takes place in intrinsic quantities only, 
that the intrinsic (mass) anisotropy ratio 7 of our PLD 
films is fully consistent with the data of the literature, 
and that the pinning constant and the critical current 
density are dominated by directional pinning in almost 
the entire angular range. The overall pinning in dc is a 
superposition of strong core pinning and Mott-insulator 
dynamic effects, and we have estimated the correspond- 
ing angular ranges. We believe that our results elucidate 
the vortex physics in presence of elongated defects. 
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FIG. 1. Microwave complex resistivity data. Full circles: real 
part; open circles: imaginary part, (a): field-sweep at 8 = 0°. 
(b): angular rotations at poHi = 0.4 T and (c) = 0.6 T. 
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FIG. 2. Flux-flow resistivity derived from Ap. (a): angular 
dependence, sample II; (b): scaling of the data over the curve 
Pff(H,9 = 0°) with the rescaled field 77//(f5); (c): experi- 
mental f(9) (symbols) and comparison with the theoretical 
3D expression, Eq.(l) (continuous line); (d) scaling for the 
field-sweeps at various angles in sample I. The log scale em- 
phasizes low-fields scaling. 
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FIG. 3. (a) Failure of the BGL scaling 1 " of k p (Hi,6) 
(open symbols) over the curve k p (H,8 = 0°) (full circles), 
(b): angular dependence of the pinning constant k p , same 
fields and symbols as in Fig. 2b. (c) Comparison between 
Jc, s hort(0.6 T,6) (see Eq.(2)) at T = 80 K and J c (l T, 6) 
at 77 K. Red arrows demonstrate the Mott-insulator effect. 
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FIG. 4. Main panel: the dynamic parameter vs. (crossed 
diamonds) shows three regions: /, where dynamic effects 
(Mott-insulator vortex phase) increase a M ; 77, where core pin- 
ning (segmented flux lines) dominates; 777, where a-b plane 
pinning is important. The green line is the empirical fit with 
Eq.(4). Inset: in region 7 the increase in depends exponen- 
tially on the component of the field along the nanocolumns. 
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